Dispersal Patterns and Influence of Air Travel During the Global Expansion of SARS-CoV-2 Variants of Concern

The Alpha, Beta and Gamma SARS-CoV-2 Variants of Concern (VOCs) co-circulated globally during 2020-21, fueling waves of infections. They were displaced by Delta during a third wave worldwide in 2021, in turn displaced by Omicron in late 2021. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of VOCs worldwide. We find that source-sink dynamics varied substantially by VOC, and identify countries that acted as global and regional hubs of dissemination. We demonstrate a declining role of presumed origin countries of VOCs to their global dispersal, estimating that India contributed <15% of Delta exports and South Africa <1-2% of Omicron dispersal. We estimate that >80 countries had received introductions of Omicron within 100 days of emergence, associated with accelerating passenger air travel and higher transmissibility. Our study highlights the rapid dispersal of highly transmissible variants with implications for genomic surveillance along the hierarchical airline network. Graphical Data analysis clarifies that dispersal of SARS-CoV-2 variants from their sites of initial detection was related to the amount of global air travel at the time of the variant's emergence, and that travel volume through "hub” sites distinct from the site of emergence was a key driver of variant spread.

Managing and assembling population-scale data streams, tools and workflows to plan for future pandemics within the INFORM-Africa Consortium

The article focuses on the establishment of the INFORM-Africa Research Hub, which aims to use big data from South Africa and Nigeria to improve pandemic preparedness. Topics include the investigation of SARS-CoV-2 transmission pathways, the impact of movement-based restrictions on COVID-19 risk factors, and the interplay between SARS-CoV-2 and HIV in Africa.

SARS-CoV-2 spike protein diversity at an intra-host level, among SARS-CoV-2 infected individuals in South Africa, 2020 to 2022.

Intra-host diversity studies are used to characterise the mutational heterogeneity of SARS-CoV-2 infections in order to understand the impact of virus-host adaptations. This study investigated the frequency and diversity of the spike (S) protein mutations within SARS-CoV-2 infected South African individuals. The study included SARS-CoV-2 respiratory samples, from individuals of all ages, received at the National Health Laboratory Service at Charlotte Maxeke Johannesburg Academic hospital, Gauteng, South Africa, from June 2020 to May 2022. Single nucleotide polymorphism (SNP) assays and whole genome sequencing were performed on a random selection of SARS-CoV-2 positive samples. The allele frequency (AF) was determined using TaqMan Genotyper software for SNP PCR analysis and galaxy.eu for analysis of FASTQ reads from sequencing. The SNP assays identified 5.3% (50/948) of Delta cases with heterogeneity at delY144 (4%; 2/50), E484Q (6%; 3/50), N501Y (2%; 1/50) and P681H (88%; 44/50), however only heterogeneity for E484Q and delY144 were confirmed by sequencing. From sequencing we identified 9% (210/2381) of cases with Beta, Delta, Omicron BA.1, BA.2.15, and BA.4 lineages that had heterogeneity in the S protein. Heterogeneity was primarily identified at positions 19 (1.4%) with T19IR (AF 0.2-0.7), 371 (92.3%) with S371FP (AF 0.1-1.0), and 484 (1.9%) with E484AK (0.2-0.7), E484AQ (AF 0.4-0.5) and E484KQ (AF 0.1-0.4). Mutations at heterozygous amino acid positions 19, 371 and 484 are known antibody escape mutations, however the impact of the combination of multiple substitutions identified at the same position is unknown. Therefore, we hypothesise that intra-host SARS-CoV-2 quasispecies with heterogeneity in the S protein facilitate competitive advantage of variants that can completely/partially evade host's natural and vaccine-induced immune responses.

Molecular Epidemiology of SARS-CoV-2 during Five COVID-19 Waves and the Significance of Low-Frequency Lineages.

SARS-CoV-2 lineages and variants of concern (VOC) have gained more efficient transmission and immune evasion properties with time. We describe the circulation of VOCs in South Africa and the potential role of low-frequency lineages on the emergence of future lineages. Whole genome sequencing was performed on SARS-CoV-2 samples from South Africa. Sequences were analysed with Nextstrain pangolin tools and Stanford University Coronavirus Antiviral & Resistance Database. In 2020, 24 lineages were detected, with B.1 (3%; 8/278), B.1.1 (16%; 45/278), B.1.1.348 (3%; 8/278), B.1.1.52 (5%; 13/278), C.1 (13%; 37/278) and C.2 (2%; 6/278) circulating during the first wave. Beta emerged late in 2020, dominating the second wave of infection. B.1 and B.1.1 continued to circulate at low frequencies in 2021 and B.1.1 re-emerged in 2022. Beta was outcompeted by Delta in 2021, which was thereafter outcompeted by Omicron sub-lineages during the 4th and 5th waves in 2022. Several significant mutations identified in VOCs were also detected in low-frequency lineages, including S68F (E protein); I82T (M protein); P13L, R203K and G204R/K (N protein); R126S (ORF3a); P323L (RdRp); and N501Y, E484K, D614G, H655Y and N679K (S protein). Low-frequency variants, together with VOCs circulating, may lead to convergence and the emergence of future lineages that may increase transmissibility, infectivity and escape vaccine-induced or natural host immunity.

Evaluation of miniaturized Illumina DNA preparation protocols for SARS-CoV-2 whole genome sequencing.

The global pandemic caused by SARS-CoV-2 has increased the demand for scalable sequencing and diagnostic methods, especially for genomic surveillance. Although next-generation sequencing has enabled large-scale genomic surveillance, the ability to sequence SARS-CoV-2 in some settings has been limited by the cost of sequencing kits and the time-consuming preparations of sequencing libraries. We compared the sequencing outcomes, cost and turn-around times obtained using the standard Illumina DNA Prep kit protocol to three modified protocols with fewer clean-up steps and different reagent volumes (full volume, half volume, one-tenth volume). We processed a single run of 47 samples under each protocol and compared the yield and mean sequence coverage. The sequencing success rate and quality for the four different reactions were as follows: the full reaction was 98.2%, the one-tenth reaction was 98.0%, the full rapid reaction was 97.5% and the half-reaction, was 97.1%. As a result, uniformity of sequence quality indicated that libraries were not affected by the change in protocol. The cost of sequencing was reduced approximately seven-fold and the time taken to prepare the library was reduced from 6.5 hours to 3 hours. The sequencing results obtained using the miniaturised volumes showed comparability to the results obtained using full volumes as described by the manufacturer. The adaptation of the protocol represents a lower-cost, streamlined approach for SARS-CoV-2 sequencing, which can be used to produce genomic data quickly and more affordably, especially in resource-constrained settings.

Genomic epidemiology of SARS-CoV-2 during the first four waves in Mozambique

Mozambique reported the first case of coronavirus disease 2019 (COVID-19) in March 2020 and it has since spread to all provinces in the country. To investigate the introductions and spread of SARS-CoV-2 in Mozambique, 1 142 whole genome sequences sampled within Mozambique were phylogenetically analyzed against a globally representative set, reflecting the first 25 months of the epidemic. The epidemic in the country was marked by four waves of infection, the first associated with B.1 ancestral lineages, while the Beta, Delta, and Omicron Variants of Concern (VOCs) were responsible for most infections and deaths during the second, third, and fourth waves. Large-scale viral exchanges occurred during the latter three waves and were largely attributed to southern African origins. Not only did the country remain vulnerable to the introductions of new variants but these variants continued to evolve within the borders of the country. Due to the Mozambican health system already under constraint, and paucity of data in Mozambique, there is a need to continue to strengthen and support genomic surveillance in the country as VOCs and Variants of interests (VOIs) are often reported from the southern African region.

Molecular Epidemiology and Diversity of SARS-CoV-2 in Ethiopia, 2020-2022.

Ethiopia is the second most populous country in Africa and the sixth most affected by COVID-19 on the continent. Despite having experienced five infection waves, >499,000 cases, and ~7500 COVID-19-related deaths as of January 2023, there is still no detailed genomic epidemiological report on the introduction and spread of SARS-CoV-2 in Ethiopia. In this study, we reconstructed and elucidated the COVID-19 epidemic dynamics. Specifically, we investigated the introduction, local transmission, ongoing evolution, and spread of SARS-CoV-2 during the first four infection waves using 353 high-quality near-whole genomes sampled in Ethiopia. Our results show that whereas viral introductions seeded the first wave, subsequent waves were seeded by local transmission. The B.1.480 lineage emerged in the first wave and notably remained in circulation even after the emergence of the Alpha variant. The B.1.480 was outcompeted by the Delta variant. Notably, Ethiopia's lack of local sequencing capacity was further limited by sporadic, uneven, and insufficient sampling that limited the incorporation of genomic epidemiology in the epidemic public health response in Ethiopia. These results highlight Ethiopia's role in SARS-CoV-2 dissemination and the urgent need for balanced, near-real-time genomic sequencing.

Subregional origins of emerging SARS-CoV-2 variants during the second pandemic wave in Côte d'Ivoire.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with increased transmissibility, virulence and immune escape abilities have heavily altered the COVID-19 pandemic's course. Deciphering local and global transmission patterns of those variants is thus key in building a profound understanding of the virus' spread around the globe. In the present study, we investigate SARS-CoV-2 variant epidemiology in Côte d'Ivoire, Western sub-Saharan Africa. We therefore generated 234 full SARS-CoV-2 genomes stemming from Central and Northern Côte d'Ivoire. Covering the first and second pandemic wave the country had been facing, we identified 20 viral lineages and showed that in Côte d'Ivoire the second pandemic wave in 2021 was driven by the spread of the Alpha (B.1.1.7) and Eta (B.1.525) variant. Our analyses are consistent with a limited number of international introductions of Alpha and Eta into Côte d'Ivoire, and those introduction events mostly stemmed from within the West African subregion. This suggests that subregional travel to Côte d'Ivoire had more impact on local pandemic waves than direct intercontinental travel.

Persistent SARS-CoV-2 infection with accumulation of mutations in a patient with poorly controlled HIV infection.

A 22-year-old female with uncontrolled advanced HIV infection was persistently infected with SARS-CoV-2 beta variant for 9 months, the virus accumulating >20 additional mutations. Antiretroviral therapy suppressed HIV and cleared SARS-CoV-2 within 6-9 weeks. Increased vigilance is warranted to benefit affected individuals and prevent the emergence of novel SARS-CoV-2 variants.

SARS-CoV-2 Genetic Diversity and Lineage Dynamics in Egypt during the First 18 Months of the Pandemic.

COVID-19 was first diagnosed in Egypt on 14 February 2020. By the end of November 2021, over 333,840 cases and 18,832 deaths had been reported. As part of the national genomic surveillance, 1027 SARS-CoV-2 near whole-genomes were generated and published by the end of July 2021. Here we describe the genomic epidemiology of SARS-CoV-2 in Egypt over this period using a subset of 976 high-quality Egyptian genomes analyzed together with a representative set of global sequences within a phylogenetic framework. A single lineage, C.36, introduced early in the pandemic was responsible for most of the cases in Egypt. Furthermore, to remain dominant in the face of mounting immunity from previous infections and vaccinations, this lineage acquired several mutations known to confer an adaptive advantage. These results highlight the value of continuous genomic surveillance in regions where VOCs are not predominant and the need for enforcement of public health measures to prevent expansion of the existing lineages.

Genomic epidemiology of the SARS-CoV-2 epidemic in Brazil.

The high numbers of COVID-19 cases and deaths in Brazil have made Latin America an epicentre of the pandemic. SARS-CoV-2 established sustained transmission in Brazil early in the pandemic, but important gaps remain in our understanding of virus transmission dynamics at a national scale. We use 17,135 near-complete genomes sampled from 27 Brazilian states and bordering country Paraguay. From March to November 2020, we detected co-circulation of multiple viral lineages that were linked to multiple importations (predominantly from Europe). After November 2020, we detected large, local transmission clusters within the country. In the absence of effective restriction measures, the epidemic progressed, and in January 2021 there was emergence and onward spread, both within and abroad, of variants of concern and variants under monitoring, including Gamma (P.1) and Zeta (P.2). We also characterized a genomic overview of the epidemic in Paraguay and detected evidence of importation of SARS-CoV-2 ancestor lineages and variants of concern from Brazil. Our findings show that genomic surveillance in Brazil enabled assessment of the real-time spread of emerging SARS-CoV-2 variants.

Emergence of SARS-CoV-2 Omicron lineages BA.4 and BA.5 in South Africa.

Three lineages (BA.1, BA.2 and BA.3) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern predominantly drove South Africa's fourth Coronavirus Disease 2019 (COVID-19) wave. We have now identified two new lineages, BA.4 and BA.5, responsible for a fifth wave of infections. The spike proteins of BA.4 and BA.5 are identical, and similar to BA.2 except for the addition of 69-70 deletion (present in the Alpha variant and the BA.1 lineage), L452R (present in the Delta variant), F486V and the wild-type amino acid at Q493. The two lineages differ only outside of the spike region. The 69-70 deletion in spike allows these lineages to be identified by the proxy marker of S-gene target failure, on the background of variants not possessing this feature. BA.4 and BA.5 have rapidly replaced BA.2, reaching more than 50% of sequenced cases in South Africa by the first week of April 2022. Using a multinomial logistic regression model, we estimated growth advantages for BA.4 and BA.5 of 0.08 (95% confidence interval (CI): 0.08-0.09) and 0.10 (95% CI: 0.09-0.11) per day, respectively, over BA.2 in South Africa. The continued discovery of genetically diverse Omicron lineages points to the hypothesis that a discrete reservoir, such as human chronic infections and/or animal hosts, is potentially contributing to further evolution and dispersal of the virus.

Identification of SARS-CoV-2 Omicron variant using spike gene target failure and genotyping assays, Gauteng, South Africa, 2021.

The circulation of Omicron BA.1 led to the rapid increase in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cases in South Africa in November 2021, which warranted the use of more rapid detection methods. We, therefore, assessed the ability to detect Omicron BA.1 using genotyping assays to identify specific mutations in SARS-CoV-2 positive samples, Gauteng province, South Africa. The TaqPath™ COVID-19 real-time polymerase chain reaction assay was performed on all samples selected to identify spike gene target failure (SGTF). SARS-CoV-2 genotyping assays were used for the detection of del69/70 and K417N mutation. Whole-genome sequencing was performed on a subset of genotyped samples to confirm these findings. Of the positive samples received, 11.0% (175/1589) were randomly selected to assess if SGTF and genotyping assays, that detect del69/70 and K417N mutations, could identify Omicron BA.1. We identified SGTF in 98.9% (173/175) of samples, of which 88.0% (154/175) had both the del69/70 and K417N mutation. The genotyped samples (45.7%; 80/175) that were sequenced confirmed Omicron BA.1 (97.5%; 78/80). Our data show that genotyping for the detection of the del69/70 and K417N coupled with SGTF is efficient to exclude Alpha and Beta variants and rapidly detect Omicron BA.1. However, we still require assays for the detection of unique mutations that will allow for the differentiation between other Omicron sublineages. Therefore, the use of genotyping assays to detect new dominant or emerging lineages of SARS-CoV-2 will be beneficial in limited-resource settings.

Comparison of SARS-CoV-2 sequencing using the ONT GridION and the Illumina MiSeq.

BACKGROUND: Over 4 million SARS-CoV-2 genomes have been sequenced globally in the past 2 years. This has been crucial in elucidating transmission chains within communities, the development of new diagnostic methods, vaccines, and antivirals. Although several sequencing technologies have been employed, Illumina and Oxford Nanopore remain the two most commonly used platforms. The sequence quality between these two platforms warrants a comparison of the genomes produced by the two technologies. Here, we compared the SARS-CoV-2 consensus genomes obtained from the Oxford Nanopore Technology GridION and the Illumina MiSeq for 28 sequencing runs. RESULTS: Our results show that the MiSeq had a significantly higher number of consensus genomes classified by Nextclade as good and mediocre compared to the GridION. The MiSeq also had a significantly higher genome coverage and mutation counts than the GridION. CONCLUSION: Due to the low genome coverage, high number of indels, and sensitivity to SARS-CoV-2 viral load noted with the GridION when compared to MiSeq, we can conclude that the MiSeq is more favourable for SARS-CoV-2 genomic surveillance, as successful genomic surveillance is dependent on high quality, near-whole consensus genomes.

Emergence and phenotypic characterization of the global SARS-CoV-2 C.1.2 lineage.

Global genomic surveillance of SARS-CoV-2 has identified variants associated with increased transmissibility, neutralization resistance and disease severity. Here we report the emergence of the PANGO lineage C.1.2, detected at low prevalence in South Africa and eleven other countries. The initial C.1.2 detection is associated with a high substitution rate, and includes changes within the spike protein that have been associated with increased transmissibility or reduced neutralization sensitivity in SARS-CoV-2 variants of concern or variants of interest. Like Beta and Delta, C.1.2 shows significantly reduced neutralization sensitivity to plasma from vaccinees and individuals infected with the ancestral D614G virus. In contrast, convalescent donors infected with either Beta or Delta show high plasma neutralization against C.1.2. These functional data suggest that vaccine efficacy against C.1.2 will be equivalent to Beta and Delta, and that prior infection with either Beta or Delta will likely offer protection against C.1.2.

Replacement of the Gamma by the Delta variant in Brazil: Impact of lineage displacement on the ongoing pandemic.

The coronavirus disease 2019 (COVID-19) epidemic in Brazil was driven mainly by the spread of Gamma (P.1), a locally emerged variant of concern (VOC) that was first detected in early January 2021. This variant was estimated to be responsible for more than 96 per cent of cases reported between January and June 2021, being associated with increased transmissibility and disease severity, a reduction in neutralization antibodies and effectiveness of treatments or vaccines, and diagnostic detection failure. Here we show that, following several importations predominantly from the USA, the Delta variant rapidly replaced Gamma after July 2021. However, in contrast to what was seen in other countries, the rapid spread of Delta did not lead to a large increase in the number of cases and deaths reported in Brazil. We suggest that this was likely due to the relatively successful early vaccination campaign coupled with natural immunity acquired following prior infection with Gamma. Our data reinforce reports of the increased transmissibility of the Delta variant and, considering the increasing concern due to the recently identified Omicron variant, argues for the necessity to strengthen genomic monitoring on a national level to quickly detect the emergence and spread of other VOCs that might threaten global health.

Rapid Replacement of SARS-CoV-2 Variants by Delta and Subsequent Arrival of Omicron, Uganda, 2021.

Genomic surveillance in Uganda showed rapid replacement of severe acute respiratory syndrome coronavirus 2 over time by variants, dominated by Delta. However, detection of the more transmissible Omicron variant among travelers and increasing community transmission highlight the need for near-real-time genomic surveillance and adherence to infection control measures to prevent future pandemic waves.

Selection Analysis Identifies Clusters of Unusual Mutational Changes in Omicron Lineage BA.1 That Likely Impact Spike Function.

Among the 30 nonsynonymous nucleotide substitutions in the Omicron S-gene are 13 that have only rarely been seen in other SARS-CoV-2 sequences. These mutations cluster within three functionally important regions of the S-gene at sites that will likely impact (1) interactions between subunits of the Spike trimer and the predisposition of subunits to shift from down to up configurations, (2) interactions of Spike with ACE2 receptors, and (3) the priming of Spike for membrane fusion. We show here that, based on both the rarity of these 13 mutations in intrapatient sequencing reads and patterns of selection at the codon sites where the mutations occur in SARS-CoV-2 and related sarbecoviruses, prior to the emergence of Omicron the mutations would have been predicted to decrease the fitness of any virus within which they occurred. We further propose that the mutations in each of the three clusters therefore cooperatively interact to both mitigate their individual fitness costs, and, in combination with other mutations, adaptively alter the function of Spike. Given the evident epidemic growth advantages of Omicron overall previously known SARS-CoV-2 lineages, it is crucial to determine both how such complex and highly adaptive mutation constellations were assembled within the Omicron S-gene, and why, despite unprecedented global genomic surveillance efforts, the early stages of this assembly process went completely undetected.